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By solving the covariant relativistic Schrodinger equations for a pair of heavy quarks, we obtained 
the wave functions for the ground and excited quarkonium states at finite temperature. In com- 
parison with the non-relativistic calculation, the J/ip dissociation temperature determined by the 
infinity size or zero binding energy of the system increases 7% — 13%, when the central potential 
varies between the free energy and internal energy. 
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It is widely accepted that there exists a quantum chromodynamics (QCD) phase transition from hadron gas to a 
new state of matter, the quark-gluon plasma (QGP) at finite temperature and baryon density. To probe the realization 
of such a phase transition in relativistic heavy ion collisions, some signatures of the new state of matter have been 
discussed for decades [l[ , and among which the quarkonium suppression is considered as a smoking gun of the formation 
of QGP0. 

A quarkonium is a deeply bound state of a pair of heavy quarks, its dissociation temperature To in a hot medium 
should be higher than the critical temperature T c for the deconfinement of light hadrons. Since a quarkonium 
is so heavy, one normally uses non-relativistic Schrodinger equation to describe its dynamical evolution at finite 
temperature. The free energy F between a pair of heavy quarks can be extracted from the lattice QCD simulations 0, 
H| ■ Taking F and the internal energy U as the two limits of the heavy quark potential V, which correspond respectively 
to a slow and a rapid quarkonium dissociation in the hot medium, the J / ip dissociation temperature Tr> determined 
by the infinity size or the zero binding energy of the system is in between 1.2T C and 2T C [1,Q. 

A nature question we ask ourselves is the relativistic correction to the dynamical evolution of a quarkonium in the 
hot medium. The correction to a bottonium is expected to be neglected safely, but for a lighter quarkonium like 
J ftp, the correction might be remarkable. The two-body Dirac equation (TBDE) f6|-FTll| of constrained dynamics was 
successfully applied to the relativistic description of light meson spectra fl2T - [l7| in vacuum. In this Letter we take the 
TBDE to calculate the charmonium wave functions at finite temperature and see the relativistic correction to the 
dissociation temperature. 

We can qualitatively estimate the relativistic effect on the quarkonium dissociation before a strict calculation. 
Neglecting the quark spin, the relative part of the Hamiltonian for a pair of heavy quarks can be approximately 
written as a non-relativistic form 



H = vV+P 2 V(r) ~ |- + V eff 
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with an effective potential 
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where /i = toq/2 is the reduced mass with mg being the heavy quark mass. Since the relativistic correction leads to 
a deeper potential well, V e jf < V, the quarkonium becomes a more deeply bound state and the temperature needed 
to dissociate the quarkonium should be higher. 

We now calculate the wave functions for a pair of heavy quarks in the frame of TBDE. Taking Pauli reduction and 
scale transformationll8j. the Dirac equation can be expressed as a covariant relativistic Schrodinger equation for a 
four-component spinor[lfJ. Explicitly, the radial motion relative to the center of mass is controlled by the following 
four equations [13(, 
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for the spin triplet uf,u^ and u 1 with quantum numbers n 2s+1 Lj — n 3 LL, n 3 Li + i and n 3 Li_i, and 
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for the spin singlet uo with quantum numbers n 1 Li, where b 2 — (m^ — 4?tjq)/4 is the energy eigenvalue in the meson 
rest frame with m m being the meson mass, n is the principal quantum number, and L, s and J are the orbital, spin 
and total angular momentum numbers. Note that the components with J = L + 1 and with J = L — 1 are 
coupled to each other. Following the notations in (l3j . we have separated the central potential into two parts, 



V(r) = A(r)+B{r), 



(5) 



and the abbreviations for the Darwin, spin-spin, spin-orbit and tensor terms introduced in the dynamical equations 
are defined as [HI 
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and m,, 



m 2 /m r , 



= (m 2 m - 2m 2 Q )/(2m m ),F' = dF/dr,G' = dG/dr and K' = dK/dr. 
We now focus on the charmonium states. The J/tp state includes two components l 3 ^! and l 3 Di determined 
simultaneously by the two coupled equations of ([3]) with energy eigenvalue b 2 = {m 2 ^ — 4m 2 )/4. Similarly, the 

two components 2 3 Si and 2 3 Z?i for the state ip' are controlled by the two coupled equations but with eigenvalue 
b 2 = (mL, — 4m 2 )/4. The xo state l 3 Po is described by the last equation of ([3]) with disappeared component , and 
the xi state l 3 Pi is characterized by the independent equation (|4} for u\. 

In vacuum, the potential between two quarks is usually taken as the Cornell form, including a Coulomb-like part 
which dominates the wave functions around r = and a linear part which leads to the quark confinement, 
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The three parameters in the model, namely the charm quark mass m c in the Schrodingcr equation and the two 
coupling constants a and a in the potential, can be fixed by fitting the charmonium masses in vacuum. By taking 
m c = 1.422 GcV, a — 0.492 and a = 0.186 GcV 2 , we obtain the charmonium masses mj/^ = 3.113 GeV, = 3.692 
GeV, m Xo = 3.404 GeV and m Xl = 3.504 GeV, which are very close to the experimental values mji^ = 3.097 GeV, 
= 3.686 GeV, m Xo = 3.415 GeV and m Xl = 3.511 GcV. 

The charmonium dissociation temperature in a hot medium can be determined by solving the corresponding dy- 
namical equation for the cc system with potential V between the two heavy quarks at finite temperature. The 
potential depends on the dissociation process in the medium. At the moment, we know only its two limits. For a 
rapid dissociation where there is no heat exchange between the heavy quarks and the medium, the potential is just 
the internal energy U, while for a slow dissociation, there is enough time for the heavy quarks to exchange heat with 
the medium, the free energy F which can be extracted from the lattice calculations is taken as the potential From 
the thermodynamic relation F = U — TS where S is the entropy, the potential well is deeper for V = U and therefore 
the dissociation temperature of charmonium states with potential V = U is higher than that with V = F. In the 
literatures, a number of effective potentials in between F and U have been used to evaluate the charmonium evolu- 
tion in QCD medium!, i, [l|l[2(|. In non-relativistic case, H.Satz and his collaborators 0, Q solved the Schrodinger 
equation and found the dissociation temperatures Td/T c = 2.1 and 1.26 for V — U and V = F respectively, where 
T c = 165 MeV[H, [f| is the critical temperature for the deconfinement. 

Considering the Debye screening at finite temperature, the potential V = F = A + B can be written asQ Q 
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where T is the Gamma function, K is the modified Bessel function of the second kind, and the temperature dependent 
parameter /x(T), namely the screening mass or the inverse screening radius, can be extracted from fitting the lattice 
simulated free energy 3- From the known free energy F, one can then obtain the other limit of the potential, 
V = U = F + TS by taking S = -OFldT. 

We use the inverse power method [21J to solve the Schrodingcr equations and (@| and obtain the radial wave 
functions 
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for the charmonium states J/ip, ip' and Xc at finite temperature. 




FIG. 1: The S'-wave (solid lines) and D-wave (dashed lines) functions for the J/ip meson at three temperatures in the limit 
of quark potential V = F. From the top down the temperature T is zero, critical temperature T c and relativistic dissociation 
temperature Td- 



The S— and D— wave functions for the meson J/-0 at different temperature are shown in FigfT] in the limit of 
V = F. The wave functions in vacuum are very similar to the results obtained in (l3| . and the component l 3 5i 
dominates the state at any temperature. With increasing temperature, both the S- and D-wave functions expand 
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continuously. At the critical temperature of deconfinement T c , while the peak values of the wave functions drop down 
a little, their behavior is similar to the one in vacuum. This means that the J/ ip meson is still a bound state at the 
deconfinement phase transition where the light mesons start to melt in the hot medium, and therefore the observed 
J ftps in the final state of heavy ion collisions can signal the QGP formation in the early stage. However, with further 
increasing temperature, the wave functions expand rapidly around T c i = 1.35T C . By calculating the average distance 
between the c and c, 



/ drr 3 \4>{r,T)\ : 
Jdrr 2 \ij){r,T)\' 



(10) 



we found that at Td the average size becomes infinite, see Fig|2j which indicates the J ftp dissociation. Therefore, Td 
is called the relativistic dissociation temperature. 
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FIG. 2: The average size of the J /if) meson in the quark-gluon plasma in the limit of potential V — F. The temperature and 
average size are scaled by their corresponding values at the critical point of the deconfinement. 



The other quantity which can be used to characterize the quarkonium dissociation is the heavy quark binding energy. 
For the Dirac equation or cquivalently the Schrodingcr equations ([3]) and (j4]), the usual binding energy defined as 
e = V^(oo) + 2m c — m m is no longer valid to describe the quarkonium dissociation. In the limit of r — > oo, the four 
equations for the spin singlet and triplet degenerate and the asymptotic equation is simplifies as 

( _ ^2 +2m w V(oo) + V 2 (oo)^ u = b 2 u (11) 

with V(oo) = -B(oo), and the energy for the scattering state is 

2m w V(oo) + V 2 (oo) = b 2 , (12) 
which leads to the binding energy for the cc bound state 

e(T) = V(oo,T) + V^ 2 (oo,T) + 4m2-m m . (13) 

FigJ3] shows the temperature dependence of the binding energy in the QGP phase and in the limit of potential 
V = F. It drops down monotonously with increasing temperature and reaches zero at Td- The result is consistent 
with the calculation of the charmonium average size. The infinite size and zero binding energy of the cc system 
define the unique dissociation temperature. In comparison with the non-relativistic calculation, the J / ip dissociation 
temperature increases from 1.26T C to 1.35T C , the relativistic correction is 7%. We also calculated the J/tp wave 
functions at finite temperature in the other limit of potential V = U and found that the dissociation temperature Td 
goes up from the non-relativistic value 2.1T C to 2.38T C . The relativistic correction is 13%. 

The wave function for the excited state xo is shown in FigH]m the limit of potential V = F. Considering the fact 
that the radial function ^ is divergent at the origin[l3|, we plot u = directly from the Schrodinger equation 
In comparison with the ground state J/ip, Xo wave function distributes in a wider region, the average size is larger 
and the binding energy is smaller. The corresponding dissociation temperature defined by (r)(Td) = oo and e(Td) = 
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FIG. 3: The binding energy of the J /if) meson in the quark-gluon plasma and in the limit of potential V = F. The temperature 
and binding energy are scaled by their corresponding values at the critical point of the deconfinement. 
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FIG. 4: The wave function u(r) = r^(r) for the xo meson at zero temperature and critical temperature of deconfinement T c 
in the limit of quark potential V — F. 



is around the critical temperature T c and a little bit higher than the non-relativistic value. It is easy to understand 
that the relativistic correction for the excited states should be smaller than the one for the ground state. 

In summary, we calculated the charmonium wave functions at finite temperature by solving the covariant relativistic 
Schrodinger equations for the cc spin singlet and triplet states. The relativistic effect makes the quark potential well 
more deep, and charmonia can survive in a more hot medium. By considering the two limits of the central potential, 
the relativistic correction to the J/ip dissociation temperature in the QGP phase is in between 7% and 13%. 
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